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ABSTRACT 


A unified method is presented for deriving the influence 
functions of moving singularities which determine the field 
quantities in aerodynamics and aeroacoustics . The moving 
singularities comprise volume and surface distributions having 
arbitrary orientations in space and to the trajectory. Hence 
one generally valid formula for the influence functions which 
reveal some universal relationships and remarkable properties in 
the disturbance fields. The derivations used are completely con 
sistent with the physical processes in the propagation field, 
such that the treatment renders new descriptions for some stan- 
dard concepts. The treatment is uniformly valid for subsonic 
and supersonic Mach numbers. 
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Geometrical Quantities 



Surface area or occupation surface 

Radial distance of a disturbance source of a field 
point laterally to a reference axis in the x direction 

Radial distance of a field point from the trajectory 
of a disturbance source (h Q = hv ) 

Length of a singularity element 

Distance between the source-sink elements of a 
dipole 


n. 


Thickness of a singularity element normal to 1 


P; P Q ; P v Location of the field point, the disturbance source 
; and the emission point 


r Q Radial distance between the field point and the 

momentary position of the disturbance source at 
time t Q 

r v Radiation radius or emission radius in spherically 

symmetrical disturbance waves 


r* 

v 


Effective emission radius when including the ex- 
tension factor of source elements 


v Q A control volume or occupation volume 

x, y, z Cartesian coordinates 





*S 


Spatial parameter of compatibility in the 
kinematics of the disturbance propagation 

Angle of inclination of the r Q line in relation to 

the trajectory of the singularities or the field 
point in the P Q PP V plane . 

Angle of emission in relation to the trajectory 
of the singularities or the field point in the 

P o PP v P lane 

Angle of emission in relation to the orientation 
of direction of the emission line or source-sink 
axis 

Angle of inclination of a source element to the 
x axis 


111 


/ 

V , 

u' ; v'; w* Components of the velocities of disturbance in 
the medium 

V„ Velocity of the disturbance source relative to 

the medium 

V„ Velocity of the field point relative to the 

medium 

V Total velocity of the medium in the disturbance 

field 

3 Mach number parameter jC^ 1 “ Ma s'l 

6 Delta function for determination of the allocation 

between disturbance source, emission point and 
field point 

e Angle of aberration between optical and acoustic 

position of the disturbance [ .0 — 0 3 

source 

T, y Basic function and unit disturbance quantity for 

the Green theorem in a disturbance field of 
moving disturbance sources, equation (4.7) 

k Relationship of specific heat [c p /c v ] 

£ Temporal compatibility parameter in the kinematics 

of the disturbance propagation 

X* Relationship between nonsteady and quasi-steady 

disturbance quantities 


p Mach angle 

p Local density of the medium in the disturbance 

field 

Density in a resting or undisturbed medium 

Ap Density disturbance in a disturbance field [p -p m ] 

a V’ °l» °n Correlation factors as a result of effective 

extension of source elements in emission 

a R Correlation factor as a result of effective 

1 shift in souce-sink elements in emission 

° D , o* Doppler factor in moving disturbance sources 

and a resting or moving field point 

<J> Disturbance potential in the medium 

Circular frequency in the nonsteady disturbance 
sources 
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v 



Angle of inclination of a source-sink axis to 

'■ the x axis 

^ Angle of inclination of the P Q PP V plane to the /7 

0 x-y plane 

Angle of inclination of the radial h<, line to 

the disturbance source in relation to the x-y 
plane 

^ Angle of inclination of the radial h Q line to the 

field point in relation to the x-hj, plane 


Aerodynamic and Aeroacoustic Quantities 


a 

a oo 

F. 

l 

G. . 
k 

Ma s 

m£ e 

m^ e 

P 

P 00 
Ap 

^co 

Q 


Local velocity of Sound 

Velocity of sound in a medium at rest 

2 2 

Disturbance force per unit of volume [kg/m sec ] 

Impulse transport parameter per unit of area 
[kg/m • sec^] 

Wave number [w/a^ = 2 tt/A] 

Mach number of the moving disturbance source 

[V/aJ 
s 00 

Mach number of the moving receiver field point 

[V E /aJ 

Mach number component of the moving receiver 
field point in the P-P Q -P V plane 

Static pressure in the disturbance field 

Static pressure in the medium at rest 

Pressure disturbance [p - p ] 

nominal dynamic pressure (reference quantity) j [| 

3 

Source intensity per unit of volume [kg/m • sec] 

Dimensionless density disturbance or pressure 
disturbance, equation (2.13) 

Singularity intensity of an element of volume, 
equation (A. 9) 


time 


t Q , t v momentary time of signal impact at the field point 

and retarded time of the emission of the disturbance 
signals 


IV 



w , u E Circular frequency of the 'disturbance signals at 

° ■' a resting or moving field point 

Disturbance function, equations (2. A) and (2.10) 
[1/sec] or [ 1 /m^] 


Indices 

o For designation of the momentary time and 

orientation between disturbance source and 
field point 

v For designation of the emission quantities for 

the impacting disturbance waves at the field 
point 

i’ j Index of direction for x, y, a (i;j = 1 ; 2 ; 3 ) 

y Emission positions (y= 1 for Ma^ci and y= 1 ; 2 

for Ma g >1) 


vi 
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I. Newton [1], J.B. Biot [2] and P.S. Laplace [3] were con- 
cerned with the theoretical calculation of sound propagation 
in homogeneous media at rest with the derivation in [2] and [3] 
completely in agreement with the physical process. The first 
formulation of the wave equation as well as the solution for the 
one-dimensional case was carried out by J.L. d'Alembert [A]. 

The wave equation was then applied extensively to propagation 
processes and vibration problems of various types. The solution 
of the radiation field with spherical symmetry in the case of 
sources fixed in space resulted from the contribution of S.D. 
Poisson [5]. The fundamental work carried out by A. Cauchy [ 6 ], 
H. v. HelmhoUTZ • [ 7 ] and G. Kirchhoff [ 8 ] on radiation fields with 
sources fixed in space of a general type are of considerable sig- 
nificance. Solutions for the two-dimensional or cylinder- 
symmetrical wave equations may be derived according to the method 
of T. Levi-Civita [ 9 ], H. Lamb [10] and J. Hadamard [11] from the 
case of spherical symmetry. In the standard texts of acoustics 
by Lord Rayleigh [12] and in hydrodynamics of H. Lamb [13], 
detained and fundamental considerations are found on the dis- 
turbance propagation of disturbance sources of different types 
fixed in space. 

The wave propagation of moving light sources was first 
studied by C. Doppler [14] and the space-time relationships in 
the disturbance propagation of moving singularities are termed 
Doppler kinematics after him. The strict mathematical treatment 
of the general Doppler principle was reported in the physics of 
W. Voigt [15]. The transformation indicated by Voigt contains 
two free constants with a suitable determination leading directly 
to the Lorentz transformation. The actual derivation according 
to H.A. Lorentz [16], or according to H. Poincare [ 1 7 ] was pro- 
duced later within the framework of electromagnetic wave propaga- 
tion of moving sources. After the introduction of the theory of 
relativity for moving systems, the derivation of the Lorentz 

Numbers in the margin indicate pagination in the foreign text. 
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transformation was presented in a more* simple and understandable 
manner, see the discussions by A. Einstein [18], H. Minkowski [19] 
as well as M. Planck [20]. Further discussions on this subject 
are found in the texts by A. Sommerfeld [21] as well as in the 
report on this topic by H. KUssner [22]. 
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The disturbance propagation of moving singularities repre- 
sents a fundamental subject in aerodynamics and in aeroacoustics . 
In contrast to the disturbance field in aeroacoustics with the 
disturbance quantities fundamentally dependent on time, there are 
non-stationary and stationary disturbance fields in aeordynamics . 
The latter represent a special case, since the disturbance fields 
or the field points accompany the moving reference system of sin- 
gularities and- the quantities of singularities themselves are not 
dependent on time. This manner of consideration is applied in 
stationary aerodymanics of moving bodies and is therefore identi- 
cal when a reference system of solid bodies and relative flow is 
employed as a basis for observations. The first reports on this 
subject, considering the effect the Mach number in stationary 
subsonic areodynamics were written by 0. Janzen.[23], Lord Ray- 
leigh [24], H. Glauert [25], L. Prandtl [26] and [27], Th . v. 
Karman [28], while the first fundamental reports on the area 
of stationary supersonic aerodynamics were presented by J. Ackeret 
[29], Th. v. Karman and N.B. Moore [30] and L. Prandtl [31]. 

A natural expansion of the subject then comprised the non- 
stationary aerodynamics of moving bodies. The first fundamental 
reports on this subject were written by H. KUssner [32], C. 

Possio [33] and I.E. Garrick [34] and [35]. 


Although the standard theory of acoustics was developed 
and expanded very early, the problems of aeroacoustics, especially 
the calculation of sound propagation from moving sources, was 
attacked much later. A simple observation and physical explana- 
tion of sound propagation of moving sources was reported by L. 
Prandtl [36]. Several important reports on this subject were 
written by H. HOnl [37], N. Rott [38]. H. Billing [39] and 
H.L. Oestreicher [40]. By converting the general Navier-Stokes 
equation into the standard wave equation for quadrupolar 
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radiation sources fixed in space, the sound propagation of tur- 
bulent air streams was formulated by M.H. Lighthill [41]. The 
solution was subsequently applied to the case of moving radia- 
tion. Further contributions on sound fields of moving sources 
were provided by I.E. Garrick and C.A. Watkins [42], H.L. Oestrei 
cher [43] and M.V. Lowson [44]. Some fundamental observations 
on energy density and energy flow in a sound field may be found 
in reports by N. Rott [38], M.J. Lighthill [41], D.I. Blokhint 
sev [45] and H.S. Ribner [46], 

When the standard methods for treating linearized field 
equations in aerodynamics or wave equations in aeroacoustics 
are considered, mainly integral methods are employed with the' 
application oT“integral transformation or the Green theorem 
with a suitable basic function. Formal mathematical transforma- 
tions analogous to the Prandtl-Glauert transformation are always 
included in the wave equation of moving sources or in the field 
equa]£^on with the Mach number terms in order to reconvert the 
equations to the known standard form. Therefore, no conclusions 
are obtained on the physical facts, resulting from the motion 
of the disturbance sources. Especially in aerodynamics, the 
effect of the Mach number is considered an effect of compressi- 
bility in the sense of a pushing together of the medium. In 
the aeroacoustics of moving disturbance sources, however, ex- 
tensive analyses of the kinematics of disturbance fields have 
been carried out, as in the reports on the subject by N. Rott 
[38] and I.E. Garrick [34] and [35]. Several further conclu- 
sions in this respect are provided in the reportsby H.L. 
Oestreicher [40] and M.J. Lighthill [41]. On the basis of the 
fact that the standard methods of solution were employed in 
these reports, there was no necessity for deriving the corre- 
lation functions of moving singularities in a general sense. 
Several fundamental observations on the origin of the correla- 
tion functions resulting from disturbance source motion were 
presented in two reports by A. Das [47] and [48]. 

In the present report, further relationships are defined 
from the basic equations on disturbance propagation and the 



A 
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resulting correlation functions from the movement of the dis- 
turbance sources are demonstrated physically. Through the in- 
troduction of the delta function in the disturbance elements, /12 

the spatial and temporal relationships may be separated carefully 
in the derivation of the correlation functions. The presenta- 
tion of the correlation functions has been generalized extensively 
by including surface occupations of singularitites in a random 
orientation in space in addition to the elementary volumes. The 
treatment and equations given apply equally to subsonic and 
supersonic velocities of the disturbance sources. The derived 
universal relationships and the resulting harmonic basic orders 
supply extensive information on this subject. Accordingly, 
several standard interpretations of aerodynamic processes in 
the disturbance field are provided in novel descriptions. 

Some fundamentals of acoustics or aeroacoustics are provided 
in the texts by P.M. Morse and K.U. Ingard [49] as well as by 
M.E. Goldstein [50]. 

2. Basic Equations and Kinematics of the Disturbance Fields of 
Moving Singularities 

In order to determine the disturbance dimensions of moving 
disturbance sources or those fixed in space, two conservation 
theorems are applied, specifically the continuity equation and 
the impulse equation, i.e. the general Bernoulli equation. By 
eliminating the quantity P through the isentropic relationship 
between pressure, density and velocity of sound, the complete 
fundamental equation of the disturbance field is obtained, 
appearing in the form of the wave equation. 

2.1. Wave Equation in Spatially Fixed and in Moving Disturbance 
Sources in a Medium at Rest 


The standard wave equation is applied to spatially fixed 
singularities in a medium at rest having infinite dimensions: 

(2.i) 7 2 ip -_Lii 2 = n (t) «(x - x ) , 

a 2 3t 2 
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2 

where <j> is the disturbance potentials A is the Laplace operator, 

ft a disturbance function and 6 the Dirac-delta function with / 1 3 

respect to space and time. 


In the case of moving singularities with the velocity Vg 


relative to the medium, the following wave equation results in 
the accompanying coordinate system, designated as the wave 
equation of moving sources. 


( 2 . 2 ) 


-L d£<£ 


,2 ( t) 6 fx 


Dt- 


oi 


W Vi 


with i = 1 , 2 and 3 for the three-axis directions and V. for the 
components of velocity in these directions. Furthermore, D/Dt 
represents th"S“substantial derivations. The position of the 
emission point (x^ , y v , ) of a disturbance wave with 

spherical symmetry and a radius r^ is clearly defined with 
respect to the momentary position P Q (s o , y Q , z Q ) with the delta 
function. It is demonstrated with this delta function that the 
path has a random orientation in space, i.e. to the three coor- 
dinate directions. Therefore, the operator assumes the following 
form for the substantial derivation: 


( 2 . 3 ) 


D 2 

Dt 2 



— ► 

+ V 



~y -v m 

with V = V„ + £v . , where v. = [ tj> , 4> , <J> ] representing the dis- 
o i i x y z 

turbance velocities in the medium. The disturbance function ft(t) 
in equation (2.1) and equation (2.2) may be produced as a result 
of the source, dipolar or quadrupolar effect and therefore' has 
the expression 

( 2 - 4 > a (t) = Q(t) + — — O-F.] (t) + — — [W-G. .] (t) . 

p » 1 "v p » l 3 

The local velocity of sound a in equations (2.1) and (2.2) has, 
with isentropic conditions, the relationship 
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with Ap = p - p ro as density disturbance. It follows from the pre 
diction of a small disturbance that: 

(2.6) a 2 = a 2 T, + Ap + 1 _ a 2 . 

« L ^ P oo J 00 

When the x axis of the coordinate system is positioned in 
the resulting path of the disturbance source relative to the 
medium and when the relationships according to equations (2.3) 
and (2.6) are introduced into equation (2.2), the linearized 
equation of the disturbance potential in a moving reference 
system is obtained after cancelling the members of higher order: 


(2.7a) 

72 


d £v> 


V V - fi 6 ( t -t - — ) 6 fx - x ■ Ma r v zl 

a 2 Dt 2 o V aj L y *o wa s y' y ' z J ' 

where ^|_ _ fa_ aj 2 

Dt 2 Lat v s 3xJ 


It also applies that 

(2.7b) 


( 1 _Ma S ) ( f xx + f yy + f zz - ^ Ma s fxt - ^ tftt = n(t)i[x v -x o -Ma s r v ,y, 2 ] 


with 


n(t) = n 6 (t -t ^) . 

o y a 


While equation (2.1) with the operator D/Dt in accordance 
with equation (2.3) supplies the complete potential equation of 
the disturbance field, the linearized equation (2.7) with the 
operator D Q /Dt is often employed in aerodynamics and aeroacoustics 
When another coordinate representation such as cylinder-polar 
coordinates, spherical-polar coordinates, etc. is introduced for 
the moving reference system, maintaining the x direction, only 
the expression appropriate for the Laplace operator is altered 
in equation (2.7) • 

In order to represent the equation of the disturbance field 
of moving sources in reference systems with fixed medium, the 
Galilei transformation is employed. As presented in Fig. 1, it 
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applies that: 


i (2.8) x = x - Vgt y = y; z = z; t = t . 


Equation (2.7) therefore assumes the following form: 


(2.9) v 2 


2 ? - *tt = 0 6( V V SI* 4 [*v-V Ma S ? v' *] 


While the left-hand sides of equation (2.1) for spatially- 
fixed disturbance sources and of equation (2.9) for moving dis- 
turbance sources have an identical form and are based on the pro- 
pagation of the disturbance waves in spherical symmetry as a re- 
sult of the compressibility of the medium, the right-hand side 
of the equation (2.9) includes the kinematic relationships, pro- 
viding the effect of Mach number in the solution. Therefore, 
the effect of the Mach number in the linearized theory for moving 
disturbance sources is a purely kinematic effect. 


When the relationship between the disturbance potential <P 
and the density or pressure disturbance s is linear, wave equations 
also result for the disturbance dimensions in the field from 
equation (2.7) or equation (2.9). These are 


2 I - 1- -2_ = o* 


v^s 


a 2 Dt 2 


Si* 5 ( t 


X 

o' V 6 [V * 0 ' Ma s f v' y> 2 ] 


(2.10) 


’ v2s " ’ 2 s tt = n * 5 <V V 5 p v - *0 ~ Ma s v y> *]■ 


In this case, 


si* = — n applies with a dimension of [1/rn ]. 


2.2. The Density and Pressure Disturbance in a Disturbance Field 
of Moving Disturbance Sources 


The density and pressure disturbances, produced at a field 
point at rest in the medium can be derived from the general 


7 



Bernoulli equation. The Bernoulli equation in a moving coordinate 
system is expressed as: 


( 2 . 11 ) i<£ + _ v s a - 

| 3t ^ - T + Th - 

The following relationship exists for the pressure distur- / 1 6 

bance in the medium: 

2 < 

(2 - 12) = E -I _ f a_] K-1 . 

P- P- UJ " 

Utilizing the relationship according to equation (2.11) for 
a/aoo , a common expression for density and pressure disturbance 
results from equation (2.12). 


In moving reference systems it applies that 

- _ A p _ A p _ 1 ph£ . „ 3jP . (?<?) 2 ~] 

s - — “ " " 77 lat v s 3x + 2 J 


(2.13) 


or in a resting reference system in the medium: 


(2.14) 


- Ac Ap _ 1 fdf , (7tp) 2 1 

s = — = Ip ^L3t 2 J 


2 2 2 2 

with ( V , ) = u* + v' + w' . This square term proves negligible 


in many cases. 

2.3. Kinematics of Disturbance Propagation in Moving Disturbance 
Sources 


When a disturbance source moves in a resting homogeneous 
medium of infinite dimensions with the velocity V s , the distur- 
bance propagation of spherical symmetry proceeds from the indi- 
vidual emision points. When a field point P(t Q ) is hit by a 
disturbance source at time t with the source located at the 

momentary position P (t ), the emission P (t ) of the wave on 

o o v v 

the path at a radial distance r = a (t - t ) from the field 

v °° o v 

point and at a distance x - x^ = Ma„r = V c (t - t ) back from 

v o s v o o v 
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from P (t Q ) (compare Fig. 2). The elapsed time t v of the emission 
is designated in comparison to the momentary time t as retarded 
time . 

When the relative position of the disturbance source P (t ) 

and of a field point P(t Q ) are described at time t by r Q and 0 Q , 

simple kinematic relationships for locating the emission position / 1 7 

P v (t v ) result for the impacting disturbance wave at the field 

point through r v , 9 y or x v and for determination of the emission 

time t . The above-described kinematics of the disturbance field 
v 

as presented in Fig. 2 is completely included in the following 
basic equations. 

(2.15) ' “ (1 - Mai) r 2 - 2r Ma cos r r 2 =0 

S V o S O V o 

(2.16) r cost^ - Ma_ r - r cos = O . 

V V b V o o 


The following expressions result as solution for the radia- 
tion quantities: 


(2.17) 


(2.18) cos 1^, 


(2.19) x 


(2.20) t 


vp 


i _Ma ! i 


| Ma g cos i^ o + (-1) w + 1 ^1-Ma 2 sin 2 t^j 


(1 - Ma|) cos 


vu 


= Ma g + 

Ma g cos t2+(-1) U+1 


= x o + Ma S r vu 


V 


1 - Mai sin 


2 tf 


fc o a 


vu 


It applies: u = i fori Ma g < i 

p = 1 ;2 for l- Ma g > i . 


It can be seen that all emission quantities determined in 
a dimensionless form according to equation (2.17) and (2.18) 
depend only on Ma^ and 0 Q . Therefore, universal diagrams may 
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be prepared for all field points in space, as presented in Fig. 3. 


The kinematic relationships in disturbance propagation accor- 
ding to Fig. 2 produce several general relationships between the 
quantities r Q , 9 o and r v > 0 V • While the field quantities in 
aerodynamics are generally described around the momentary source 
position P Q by means of r Q and 0 , in aeroacoustics a representa- 
tion is generally employed around the emission position by 

means of r y and 8^ . In this case, the following relationships 
apply: 


/ 1 8 


( 2 . 21 ) 


r v u (cos* v -Ma ) 


r vu sin % 


- r cos ■>> 

o o 


= r sin v' 
o o 


(- 1 ) 


u+1 


r vu (1 ~ Ma s cos = r o n-Ma| sin 2 ^) 1 ^ 


r vu (1 + Ma | - 2 Ma costf ) 1/2 = 


v u 


sin e 


= Ma c sin . 
S o 


Furthermore, the following expression is obtained from 
equation (2.19) for Ma s >1 with P = 1 and 2: 


( 2 . 22 ) 


X V2 ~ X VI = ^S 4 _ 2 Ma s ■%/ 1 ~ Ma i sin 2 '?,' 

r " r lt vi ~ t V 2 ) = ! • 


Ma | ~ 1 


The relationship always applies for the propagation of 
disturbance waves in spherical symmetry: 


( 2 . 23 ) 


(X - X u ) 2 + (y - y v )2 + (2 - 2^)2 _ a 2 (t _ 2 . Q . 


Equation (2.23) remains invariant in the transition from a 
reference system at rest to a moving one, as well as in the coor- 
dinate transformation from physical space in the Lorentz space. 

When a source element is shifted from the disturbance center 
P 0 (x 0 , y , z ) along or across the path line, a corresponding 
shift of the emission point occurs according to the relationships 
in equations (2.17) to (2.21) under the condition that the 
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disturbance signals reach the field point at time t Q . The fol- 
lowing relationships in the differential form have proven very 
useful . 


In the shift of a disturbance source along the path line, 
the following relationships result: 


3x v 


1 


3x o 

(i 

- Ma g cos 

V 

3 \ 

3x 

o 

0 



3r 

v . 


COS 
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3x o 

(1 

- Ma g cos 

V 

3*0* 

V 


sin *0* 

V 



V 3x o (1 - Ma„ cos ) 

b V 
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When a shift dh^ of the disturbance source is undertaken in 
a random direction in the space across the path line with the 
field point located at P(x, h Q , T 0 ), the relationship applies: 


( 2 . 25 ) 

with 

( 2 . 26 ) 



COS ijl 


h = r sin -\> = r sin = h . 

o o o v \> v 


In the shift of a disturbance source across the path line, 
the following expressions therefore apply: 


( 2 . 27 ) 
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( 2 . 27 ) 


3i)< 

r 

r v 3h s 


(cos - Ma g ) cos p 
( 1 - Ma s cos iJ^ ) 


The relationships according to equations (2.24) and (2.27) 
are used later for determining the correlation functions. 

3. Derivation of the Correlation Functions in Disturbance Fields 
of Moving Singularities 

In order to determine the disturbance quantities at a field 
point at time t Q , the radiation radius r v of the impacting dis- 
turbance wave is required, of course, the strength ft(t v ) of 
the disturban ce source at the retarded time t of the emission 
and all the correlation functions, resulting from the motion of 
the disturbance source. As discussed in detail in [47], the 
correlation functions are derived from three physical observa- 
tions : 

a) the effective extension of source elements in emission, 
with the disturbance signals simultaneously reaching 
the field point at time t Q . The effective extension 
results in a longitudinal or transverse direction, 
therefore producing a volume expansion. 

b) the effective shift of source-sink elements in emission 

with the disturbance signals simultaneously reaching the 

field point at time t . 

o 

c) the standard Doppler effect as a result of the altera- 
tion in wave passage duration dt Q through the field 
point in comparison to the duration dt^ of the emission 
of the disturbance signals. 

All these correlation functions may be derived directly from 
the basic equations, (2.17) to (2.27). 

The following designations are introduced for the corre- 
lation functions: 
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The Spatial Correlation Functions 


dl 

\i 

°L = dT 


dn 

_ v 

3 N - dn 


dv 

v 

dv 


correlation factor as a result of effective extension 
of the emission elements in a longitudinal direction 
coinciding with the Mach number component Ma g = 

Ma s cos Xg. 

correlation factor as a result of effective extension 

of the emission elements in a transverse direction 

coinciding with the normal components of the Mach 

number, Ma = Ma~ sin x o . 

n o s 

correlation factor as a result of effective extension 
of_ the volume elements of the disturbance sources in 
emission . 


dr v 
a R = dF 


VO 


correlation factor as a result of effective shifting 
in the source-sink elements during emission. 

(dr vo = dr v at Ma S = 0) ' 


The Temporal Correlation Functions 


= Doppler factor in medium-fixed field point 


D dt 


dt 

_* _ V 

D dt 


Ha, 


Doppler factor in moving field point. 


These correlation factors are now derived for source elements 
with random orientation in space and to the path line. 


3.1. The Spatial Correlation Functions as a Result of Effective 
Extension and Shifting of Source Elements in Emission 


The Correlation Function ° ^ 

When the disturbance effect of a moving elementary source 
path dl Q is observed in space with random slanted position x s 


/2 1 
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to the path line, a correlated emission path dl v results accor- 
ding to the kinematic relationships of section 2. This extension 
in length occurs in the direction Ma Q = Ma^ cos Xg . The process 
described is illustrated in Fig. 4. The effective emission path 
is obtained by combining the effective partial elements <J>1 , 

with disturbance signals reaching the field point P simultaneously 
at time t . It then follows 


( 3 . 1 ) 


dl 

V 

dl 

O 


16 1 
dl o 


dx dh 

dT cos *s + dl~ Sin X S * 
o o 


Through utilization of the partial derivations according to 
the relationships in section 2, the following exact equations for 
the correlation functions result. It then applies that 


( 3 . 2 ) 



1 + 


Ma cos •»> 
e ve 

T - Ma_ cos v v 
S v 


( 3 . 3 ) 


dl 

v _ 1_ 

1 - Ma* cos i>* ' 
e v 


where Ma 0 = Ma s cos Xg an< ^ = ^ a c / cos (X s - X v ) representing 
the Mach number components along dl v and dl^. 


The Correlation 


Function 


When a moving source element exhibits a random orientation 
in space or to the path line, an altered transverse extension 
dn v of the emission element is obtained in comparison to the 
actual width dn Q from the same observation as in the previous 
section. This fact is illustrated in Fig. 4. It follows from 
the kinematic relationships in this figure: 

dn 

(3.4) 0 = v = 1 

N dn 1 - Ma„ cos 

o n vn 

i 

with Ma n - Mag sin Xg, representing the Mach number component 
in the direction normal to dl Q and forming 9 Vn the emission angle 
between Na n and r v 
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The Correlation Function a v » 

When a moving element of volume of disturbance sources is 
observed with random orientation in space and to the path line, 
the effective emission volume dv^ results with disturbance signals 
reaching the field point P simultaneously at time t Q . On the 
basis of the fact that the element is not extended normally to 
the direction of motion Ma s , the alteration in volume consists 
of the longitudinal and transverse extension as already described. 
It therefore follows that 

j dv y dl y dn v 

j ( 3 . 5 ) °V dv Q dl Q * dri^ “ °L " °N 


and may apply— that 




Ma„ cos - 0 * ) 

( 3 . 6 ) 

a 

1 + e ve 1 

V 

1 - Ma_ cos 1 

S v > 


From the relation 


i 

1 - Ma cos 

n vn 


( 3 . 7 ) Ma_ cos •«> = Ma cos ^ + Ma cos ■<> 

° v = ve n vn 


it follows for the effective volume extension of a moving source 
element with random orientation in space: 

i 

1 - Ma 0 cos 

S v 

The type of origin of the correlation function is illus- 
trated in Fig. 5 for various orientations of an element of volume. 



This results in 


°v = ° L 


1 - Ma c cos •»> 

b V 


a N = 1 


( 3 . 9 ) 


°V = ° L • ° N = 


1 - Ma cos 

b v 


°v " = 


1 - Ma_ cos ■»> 
b v 


°L = 1 


for x s = o 


for : 0< x s < 7 


for. x s 


If 

2 • 
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It is noteworthy that the volume expansion and therefore 
in the present three cases has the identical value. 


The Correlation Function 

When a shift is undertaken from the amount dm Q of the source 
element at P Q in the direction 9^ to the path line, a shift re- . 
suits at the emission point P y for the source- sink components 
with signals simultaneously reaching the field point at time t Q . 
This situation is illustrated in Fig. 6. 


In a similar procedure to that in the previous section, the 
difference in emission radii of source-sink elements, following 
a path at a di-stance dm Q , can easily be determined at a field 
point P. It applies that: 


( 3 . 10 ) 


dr 

v _ 

3r 
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3r 

v . 

3h o 

dm 
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3x o ’ 

3ni 

o 

3h 

0 

3m 
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By employing the partial derivations according to the rela- 
tionships in section 2, the exact expression is obtained: 


( 3 . 11 ) 


dr 

V 

dm 


cos iV 
v 


o 1 - Ma n 'cos 

S * v 


where 0 v represents the angle of inclination between r v and dm Q 
at the emission point P v and exemplifies the following relation- 
ship . 


( 3 . 12 ) 


COS *5 
V 


cos cos Xg - Sin -0^ cos(<J> 0 - * s ) sin x g . 


When the actual source-sink arrangement at P Q is applied as 
spatially-fixed system at the emission point P v , the difference 
results there in the radiation radii of these elements to the 
field point P at 

i 

(Py) = COS . 

Ma s =° 


( 3 . 13 ) 


dr 
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The correlation function a R follows from equations (3.11) 
and (3.13) and supplies the generally valid expression: 


(3.74) 


dr /dm 
v o 


| dr /dm 
1 v o 



= 1 

dr 1 - Ma_ cos 

vo S v 


The expression o R according to equation (3.14) represents a 
universal formula and is applied for elements of volume as well 
as for surface occupation of singularities. 


3.2. The Temporal Correlation Functions Resulting from the Altera- 
tion in Wave Passage Duration at the Field Point in Comparison 
to the Emission Duration of Disturbance Signals 


In the cffsturbance propagation of moving disturbance sources, 
the wave passage duration dt Q generally differs from the emission 
duration dt^ of the disturbance signals by a field point. There- 
fore, a temporal extension effect dt /dt Q , the so-called Doppler 
effect, arises, defined in the case of medium-fixed field points 
by a D and in the moving field point by ag. While the correlation 
factor a D also determined the field dimensions in the medium, 
is employed for frequency determination of the received signals 
at the moving field points. This situation is illustrated is 
Fig. 7. The derivation of equation (2.20) for t results in 
the universal expression: 


dt o 

(3.15) 


3r /at 
o o 

a.(1-Ma|) 



3^ Q /3t 
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( 7 -Ma^) 


Ma_ cos 


+ 
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(- 1 ) 


11+ 1 




-MaJ 


sin 2 !?- 




Ma sin i* + (- 1 ) P + 1 

o O 


Ma 2 sin 2 i> 

S o 

2 -^1-Ma 2 sin 2 -^ o ' 


The temporal alteration of r Q and 


(3.16) 



(Ma_ cos ifr 
S o 


r 3i> 

_£ 2 = (Ma_ sin •»> 

a 3t S o 

a* O 


0 o is given by: 


+ Ma_ cos ) 
t. o 


+ Ma sin »> ) . 
E o 
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In this case Ma^. is the component. of the field point Mach 
number relative to the medium in the P Q PP^ plane and e Q is the 
angle between the vector Ma^ and r Q in the same plane. In the 
case of a medium-fixed field point, Ma^. = 0. 

The Correlation Function a D 

From equations (3.15) and (3.16) the result for the Doppler 
factor at a medium-fixed field point: 


(3.17) 

dt 

‘"d = St* 

O 


.1 

— 

1 + (-1) p+1 

Ma cos 

o 




1 ’ Ma s 




Ma E =0 


J 1 -Ma 2 sin 2 ^ 

* b o 


The relationships according to equation (2.21) result in: 


(3.18) 


°D = <- 1} 


u+1 


( 1 - Ma_ cos ) 
S vii 


with »■ =i 

u =1 ; 2 


for Ma s < 1 


for 


Ma s > 1 


The Correlation Function 


In a moving disturbance source and moving field point with 
a finite value of Ma E , the universal Doppler factor o* f results, 
decisive for the frequency determination of disturbance signals 
at the field point. According to equations (3.15) and (3.16) it 
applies that: 


dt 
v 

d t 


■Ma r 


1 - Ma| 


Ma„ Ma_ cos (i* + -O' ) 

b ti O O 

1 - Ma* 


(3.19) 
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where Mlg and F Q are explained in Fig. 2. It is often practi- 
cal to rewrite the Doppler factor a* in the expression of 0 and 
6" . For this purpose, the relationships according to equation 

(2.21) are employed as well as the following relationships: 

I . r o cos = [ r v cos5 v - Ma s r v C os(^ v + ?)] 

(3.20) u 

"1 • 

o o L v vj „ 

The conversion of the expression of a* according to equation 
(3.19) to the reference dimensions at the emission point supplies: 


u 


D 


(3.21) 


. Ma Ma cos(i) +•»> ) 

1 + S E v v + 

1 - Ma| 1 - MaJ 

Ma s (cos - Ma g ) + Ma E | cos - Ma g cos (^+7) 

(1 - MaJ) (1 - Ma„ cos ■»> ) 

o v 


Ma| Ma sin -O' sin( %> + -»>) 

_ a ts V V V 

(1 - Ma|) (1 - Ma cos ) 
t> S v 

with u =i for Ma s < 1 

and , ’ u =i ; 2 for Ma g > i . 


When Ma E = 0 is placed in this equation, then a* converts 
into the expression of o D according to equation (3.19). The 
general expression of the Doppler factor o D according to equa- 
tions (3.19) and (3.21) includes all simple cases of the moving 
field points or receivers, occurring in literature. A repre- 
sentation of the results according to equation (3.18) for sub- 
sonic and supersonic Mach numbers of the disturbance source- is 
carried out in Fig. 8. It is noteworthy that the Doppler factor 
o p for 0 ^ = 0 assumes the maximum value at Ma^ = 0 and always 
becomes smaller with decreasing or increasing Mach numbers of the 
disturbance source, as the figure on the right demonstrates. 

While the Doppler factor constantly remains at 1 in the case of 
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0 = tt/ 2 , a continuous decrease in the a n values is recorded 

v D 

for e v = 7T with increasing Mach number. 


3.3. Universal Representation of the Correlation Functions and 
Several Special Characteristics in Disturbance Propagation 

The spatial and temporal correlation functions, arising as 
the result of the disturbance source motion, may be compiled in 
a universal representation for random field points in space. When 
the correlation functions a^, and a D for random medium-fixed 

field points in space are plotted as a function of Ma s cose^ , 
a universal curve results, consisting of two hyperbolic branches 
valid for a Mach numbers and all field points. This is explained 
in Fig. 9 . The appropriate range of validity of the curve is 
extended by the value a= 1 for Mag = 0 with increasing Mach 
number, as can be seen from this illustration^ In the case of 
Mag > 1 , the range extends beyond the line Ma g cose^ = 1 , re- 
presenting the point of contact of the wave having spherical 
symmetry with the Mach cone. When a wave having spherical 
symmetry is observed during propagation around an emission point 
, the correlation factors a in this wave remain at the same 
value for each 0 V = const. 


When the correlation factor is combined with the radia- 
tion radius r v , a further concept is created in the form of 
the actual radiation radius according to the definition r* = 
r v / a y. It then follows according to equation (3.8): 


( 3 . 22 ) 



r (1 - Ma_ cos ) . 
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The alternate expression then obtained according to the 
relationship in equation ( 2 . 21 } 


( 3 . 23 ) 

top. 


r Q (1 - Ma| sin 2 ^)^ 2 
( 1 - Ma 2 sin 2 i> ) ^ . 
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A universal curve, i.e. an ellipse, is then obtained for all 
Mach numbers and for field points in the entire space, as is illus- 
trated in Fig. 10. The range of validity of this curve increases 
with Mach number, achieving the limit value r* = 0 for Mag > 1 

at Ma s sin 0 q = 1. An important characteristic can be seen in 
Fig. 10, as the quantities r v /r Q , 0 v and r* /r Q always remain 
constant on a radial line of the momentary source position P , 
i.e. for 0 = constant. It therefore follows that the impacting 

disturbance waves always have the same correlation functions 
at field points on this line. 


The radiation process in an evenly moving source path of 
infinite length and random inclination to the path line is illus- 
trated in Fig. 11. In this case, a field point P is observed on 
the same plane as the source path when the disturbance signals 
impact there at time t Q . The emission points of this signal can 
easily be determined according to equation (2.17) or equation 
(2.19). Since the position of the source path at time t Q is known, 
the position and shape of the emission line L may be completely 
defined from the guideline relationships. It applies that 


( 3 . 24 ) 



r 


v 

Ma r 
S v 



This equation supplies a hyperbola for Ma s < 1, a parabola 
for Ma^ = 1 and an ellipse for Mag > 1, where the field point is 
correlated as focus of these curves. For random field points 
in space, corresponding lines are obtained, the envelops of 
these base curves. The base curves are altered with respect to 
position and extension with the slanted position of the source 
path to the path line by assuming a symmetrical course around the 
slanted axis. The methods for solving the disturbance fields 
in standard area dynamics are essentially supported on this 
characteristic, by basing the solutions on the Mach number com- 
ponent Ma n = Mag sin X of the source path. As is clearly demon- 
strated in Fig. 11, this procedure is an infraction of physics 
at Ma s < 1 , since a portion of the disturbance sources situated 
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outside of the actual Mach cone is considered effective. The 
situation can easily be seen in Fig. 11. Therefore, it is clear 
that the division of the motion Mach number Ma s into the normal 
and tangential components Ma n and Ma e , as is undertaken in 
standard aerodynamics, is a completely artificial method to 
achieve a solution. The actual sweepback effect of a slanted 
source path can be represented as follows. 

In moving sources paths of infinite length, the actual 
extension of the emission elements in the longitudinal direction 
dl Q has no additional effect, since the entire length is already 
occupied by the sources. In contrast, the actual extension of 
the emission elements in the transverse direction is completely 
effective and—is included entirely by the correlation function 
a . This situation is in complete agreement with the physical 
processes in the disturbance field as shown in Fig. 4 and supplies 
the true basis for the sweepback effect. 

A further noteworthy process in disturbance propagation of 
moving source paths in the case of supersonic Mach numbers is 
demonstrated in Fig. 12. In this case, it is shown that the 
disturbance waves impacting at the field point P in each case 
at times t Q ^ and originate at the emission lines L,^ and L^. 

It is remarkable that the entire effect of the disturbance signals 
remains constant at the field point for all times t Q , since the 
radiation process with r , a and L always supplies the same 
amount for the surface integrals of the 1/r* curves. It can be 
demonstrated that the conditions remain in all field points in 
space behind the source path, specifically under the influence 
of the above-described dynamic processes and continuous expan- /3 1 

sion of the emission lines. 

When a source path of finite length moves in a longitudinal 
direction with the source elements situated along the path line 
X S = 0, as is the case in axis-symmetrical disturbance fields, 
several remarkable phenomena occur in the propagation process 
along the path line. When a field point on the path line of the 
source path with X s = 0 is observed, the effect of Mach number 
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disappears there with the effective emission radius r^*. The rea- 
son for this lies in the fact that in the appropriate 0 O = 0 or 
tt of the extension factor cry of a source element and the dimension- 
less emission radius r /r of the disturbance wave at the field 

v o 

point assume identical values as is demonstrated in Fig. 13. Al- 
though both quantities and r^ are subject to the entire 
effect of Mach number, the total effect is cancelled at r* with 
the result that r* is identical to r Q . It applies according to 
equations (2.21) and (3.13): 


r 

(3.25) — ■ = o v 

o 


1 

1 iMa s 


Therefore it follows that: 


for i> 0 = o , n . 


(3.26) r* 

V 



for = O T>d n . 

O r 


This situation supplies the physical basis for the theory 
of thin bodies, postulated in literature by other means. 


When three-dimensional singularity occupations of any type 
are located in smooth motion in space, all singularities with r* = 
constant may be combined for determining the disturbance quantities 
at a field point at time t Q . The momentary position of these sin- 
gularities in relation to the field point are described by r Q and 
0 , necessarily fulfilling the relationships according to equa- 
tion (3.23). It follows that: 


(3.27) 


( 1 - Ma 2 sin 2 i> ) 


vr 


In the case of rj = const, and Ma s = const., r Q describes 
surfaces or shells for the following form according to this 
equation : 

a rotational ellipsoid for Ma s < 1 

a vertical plane for Ma^ = 1 

a rotational hyperboloid for Ma<, > 1 , 

where the x axis forms a symmetrical axis through the field point. 
These functional characteristics supply the basis, designated in 
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literature as surface rule. Several examples of this are calcu- 
lated for three Mach numbers from equation (3.27) and illustrated 
in Fig. 14. 


In aerodynamics and aeroacoustics , problems in surrounding 
flow or radiation are often simulated with the surface occupa- 
tion of singularities. When the correlation functions are deter- 
mined for this purpose as a result of the motion of the surface 
elements, the singularity elements must be considered as elements 
of volume of thin shells, because physical line occupation or 
surface occupation of the thickness of zero cannot be realized. 
The extension effects in emission are supplied by the correla- 
tion function according to equation (3.9). According to this 
it applies tharfc: 


(3.28) 



i 


1 - Ma_ cos 

S v 


The analysis of this expression with the aid of equation (3.9) 
supplies some important new knowledge. 


While a pure surface extension of the disturbance elements at 
Xg = 0 is recorded with = cr and = 1 , an alteration in flow 
amount or source density results at Xg = tf/2 through the emission 
surface dl„ = dl , because cr = i remains. This situation is 
caused by the fact that no surface extension occurs in the ele- 
ment but only a transverse extension through = a^. In the 
case of 0 < Xg < tt/2, both effects occur together with emission. 
This situation is demonstrated in Fig. 5. According to the 
above-described physical process, the following relationships 
result for moving surface elements in a random orientation to 
the path line: 

S = a S 
v N o 

j(3 . 29) 


According to equation (3.28) it also applies that 

(3 ' 30) °v s o dF o = S v dF u • 
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Further fundamental analyses of the disturbance fields may 
be undertaken with the derived kinematic relationships and cor- 
relation functions, and many standard interpretations may be 
formulated in accordance with the facts in a new representation. 

4. Solution of the Wave Equation and Derivation of the Field 
Quantities in Moving Singularities 


The propagation processes of spatially fixed and moving 
disturbance sources in a medium at rest of infinite extension 
may be considered in the same manner as medium-fixed or moving 
reference system. The medium always serves as a carrier for 
the wave system with the disturbance centers and field points, 
and therefore the wave image, remaining identical in both systems. 

The linearized wave equation in the moving reference system has 
the general expression: 

i d 2 <# 

(4.1) 7 2 * - ^ — = n 5(< ) 5(e ) . 

a 2 Dt 2 ° ° 

“ 

In the medium-fixed reference system, the wave equation 
assumes the following form: 

« 

(4.2) v2 -f " h ^ = " S( *o } 5(c o } • 

ai 3t 2 

The delta functions in these equations include the temporal 
and spatial kinematic relationships in a universal form, since 
they apply both to spatially fixed and to moving disturbance 
sources. The designations introduced in this represent the /34 

following relationships: 


(4.3) 



r 


V 


In spatially fixed disturbance sources, Ma s = 0, while the 
motion Mach number of the singularity is ft in moving disturbance 
sources Ma s> The expression for £ Q is identical in both cases. 

As already described in section 2.1, the equations (4.2) and (4.3) 
can be converted to one another through the application of the 
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Galilei transformation. The delta functions in both equations 
remain invariant in this process. 

In order to solve the wave equation of moving disturbance 
sources, generally integral methods are applied, utilizing in- 
tegral transformations or through the application of the Green 
formula with a suitable base function. In both these methods, 
first the field equation is converted into the standard form. 

For this purpose a mathematical transformation analogous to the 
Prandtl-Glauert transformation or a Lorentz transformation is 
employed. In these transformed spaces, the equivalent radiation 
radii assume the same dimensions ,• identical to in equation 
(3.23). In standard aerodynamics, this situation is attributed 
to the effect— of compressibility. 

In the following section, the integral method in the actual 
physical space is employed for solving the wave equation of 
moving disturbance sources, specifically in such a manner that 
the correlation functions are included appropriately in the 
Green function, i.e. in the base function. 

4.1. Solution for the Wave Equation According to the Expanded 

Green Theorem 


In order to treat the wave equation in the case of moving 
singularities, the two expression according to equations (4.1) 
and (4.3) can be equally applied. For the subsequent suitable 
steps, the equation is employed here in a medium-fixed reference 
system. It states: 


( 4 . 4 ) 



n & U o ) 6 (e q ) . 
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In comparison to the wave equation of spatially fixed sin- 
gularities, two important characteristics must be underlined 
here, included in the solution. The motion of the disturbance 
source produces a shift in disturbance centers and causes a 
spatial extension effect of the emission elements with the 
singularities n. These two effects, determining the disturbance 
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potential at the field point, are always completely included in 
g and (Mag, r^ , ej . 

The disturbance potential 4 ) (t) dependent on time at the 
field point may be constructed by means of a Fourier series from 
a number of harmonic functions with the form 

— ioj t 

(4.5) 4 >(t) = *f Q e 


In order to solve the wave equation by means of the Green 
theorem, a basic function is now introduced such that it satisfies 
the wave equation for harmonic disturbance functions and includes 
the spatial correlation created at the emission center. The base 
function is expressed as 


(4.6) 


r (p) 


i_ ^v 

4tt r 
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.ikr 
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with k = m /a m , r v as radiation radius of the disturbance wave 
at the field point. Furthermore, o ^ according to equation (3.8) 
represents the extension function of the correlated emission 
element at the disturbance center. While the correlation func- 
tion a v (Ma s , 0 v ) is not dependent on r v and V , this assumes a 
fixed value corresponding to equation (3.8) for an infinitesimal 
surface element of the disturbance wave at field point P(r y , 0 V ) , 
because the actual emission element for the disturbance wave at 
the field point is completely defined. Therefore, the basic func- 
tion fulfills the Helmholtz equation in the following form: 

(V 2 + k 2 )r = y£(<) 5(e) 

o o 

with Y as a unified disturbance function. When the same proce- 
dure is applied as according to Helmholtz [ 7 ] and Kirchhoff [ 8 ] 
and the delta function <M^ 0 ) is introduced into the distur- 
bance quantities in an integral form in the application of the 
Green theorem, the solution of the disturbance potention results 
in the following form: 
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where p signifies the number of positions on the path line, ful- 
filling the condition t Q = 0 and n Q representing the surface 
normal value for F . It follows from section 2 that p = 1 for 
Mag < 1 and = 1 and 2 for Ma^ > 1 . In spatially fixed dis- 
turbance sources with Ma s =0, - 1 and, as a result, 

equation (A. 8) is converted to the solution according to Helmholtz 
[7] and Kirchhoff [8], The wave equation of moving disturbance 
sources may also be solved in an elegant manner by the introduc- 
tion of the Lorentz transformation. This subject will be treated 
in more detail elsewhere. 

The first integral in equation (4.8) includes the amount of 
the element of volume of singularities with the strength ft (t), 
while the second integral represents the amount resulting from 
the source and dipolar occupation on the surface F . 

4.2. Determination of the Disturbance Potential in Moving 
Singularities of Various Types 

The solution of the wave equation of moving disturbance 
sources with their general form according to equation (4.8) 
includes the disturbance function ft(t) in the volume integral, 
where these may represent a source, dipolar or quadrupolar sin- 
gularity according to equation (2.4). The dipolar and quadru- 
polar singularities may in turn assume various arrangements in 
space or in relation to the path line. 

When the three disturbance functions in equation (2.4) are 
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considered harmonic functions, where Q^, F q and G q form the maxi- 
mum amplitude of the disturbance dimensions, the following inte- 
gral quantities are introduced in the volume integral: 


i • § ° = k J Q o dv o = M o 

(sec 

(source) 

(4.9) V ° 



S 2 = — I F dv = D 

w v p « I o o o 

V 

o 

m 4 

sec 

(dipole) 

/ G oy ■>'<, ■ T o 

"m 5 I 

sec] 

( quadrupole) 


V 

O 


A general solution of the volume integral in equation (4.8) 
is obtained by— the application of partial integration resulting 
in a unified expression of the following form: 


(4.10) 


(-D 


u+1 


l+m 
4 IT 


. l+m 


ax* ah': 
o o 


-iu t , 
v VP 


VP 


(1 - Ma c cos ^ ) 

b V U 


(with a V p = 1/(1 " Ma s cos ^.p ) ' 

where the following combinations result for the disturbance 
functions : 

Sources: 1 = m = 0 

Dipole: l=1,m=0;l=0,m=1 

Quadrupole: 1 =2, m=0; 1 = 1, m = 1; 1=0, m = 2. 

The surface integrals in equation (4.8) may be calculated 
analytically or numerically according to the singularity occu- 
pation with the surface being divided in the latter case into 
individual elements. The determination of the induced effect 
of these surface elements at the field point may be carried out 
such that the integral expressions are first described in the 
form of contour integrals along the edges of the surface elements. 
The desired field dimensions result from this. 
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4.3. Determination of the Density and .Pressure Disturbances 
in the Disturbance Field 


The field quantities with respect to density and pressure 
disturbances can be calculated with a unified formula according 
to equation (2.14). It applies in a medium-fixed reference 
system that: 


(4.11) 


I(P) = ± 2 . « 


AP = _’L_ 




where the functional relationship of <f> = <f>(Maoj r M > 6,, > t ) 

O V V V 

is included in the derivation of <J> for t . The determination of 
_ ° 

s can be carried out in a few steps, if the appropriate delta 
function 5(£ ) is introduced in integral form in the derivation 
according to equation (4.11). It then applies that: 

CO 

ft; | {* ( V <o ) 4 <*o>} d "o + • 

“ ‘ K 0 -*° J 
After carrying out the derivation for t and after partial 
integration of the delta function, it follows in the moving re- 
ference system that: 


(4.12) 


i (P) = Ae. = A£_ = - JL 
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or 

(4.14) 
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The derivations achieved when fulfilling the kinematic re- 
lationships £ D = 0 and £ q = 0 were already derived in sections 
2 and 3. The field quantities are clearly described there and 
the resulting terms reproduced as near field and distant field 
portions. Exact relationships between instationary, quasi- 
stationary and stationary field quantities may be prepared from 
these figures. 
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In a moving source volume, the following density and pressure 
disturbances occur in the disturbance field: 

Near Field Term 


(4.15) 




(-D 


u+1 


S(t ) 
vu 


Ma s cos 


( 1 -Ma| 


sin 


^J 3/ M 


Distant Field Term 


^ ( 1 - Ma cos oS ) I 
i> v ' 

Equations_(4 . 15) and (4.16) produce the instationary dis- 
turbance quantities in the field. The disturbance quantities 
and Sp on the unit spherical surface around x Q and x y are 
illustrated in Fig. 17. When Mag<< 1, these quations convert 
into the quasi-stationary forms by replacing the expressions in 
parentheses by 1 . The relationship of the instationary to the 
quasi-stationary disturbance quantities produces 


s f * 




u+1 iw S(t ) 

V vu 


(4.16) 

and 


i 

(1 -Mai sin 2 ) 3/<2 
S O 


(near field)' 


(1 - Ma s ccs x? r 


(distant fi-eld) 


These factors are reproduced for various Mach numbers and 
all field points in space in Fig. 18. They supply some interes- 
ting conclusions with respect to the effect of Mach number and 
field point orientations with the field points located on beams 

with constants 0„ or 8 

o v 

In the temporal, constant source intensity, stationary 
fields result with accompanying reference system by setting S = 
S Q = const, in equation (4.15). It can be seen that the distant 
field portion of the field quantities disappears in this case, 
because the quantity oj v = 0 is the equation for Sp. When a 
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moving disturbance source is observed in an incompressible medium, 
it is practical to replace the quantity s^ by the coefficient of 
pressure Cp. It then applies that 


(4.18) 



2_ 
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4tt r 2 
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In the case of instationary disturbance fields, the deter- 
mination of the signal frequency at the resting or moving field 
point is carried out such that the emission and reception of a 
defined number of waves is taken into consideration. It is 
necessary that 


(4.19) 


0) dt 
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v v 
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Subsequently 

I (4.20) 
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a io 
D v 


and “e ~ ? d“v * 


The temporal extension functions ci D and a* are already known 
from equations (3.18) and (3.19). 


5. Summary 

First all the correlation functions resulting from the motion 
of singularities are derived in a generally valid form, using the 
kinematics of disturbance propagation from moving disturbance 
sources by permitting random orientation in space, i.e. to the 
path line, for the source elements. It follows that the resulting 
correlation functions have spatial and temporal extension effects 
in the emission of the disturbance signals. 

The correlation effects may be obtained in a universal 
representation for a random number of field points in space and 
for all Mach numbers in individual diagrams. All physical processes 
in the disturbance field are factually included by the the present 
manner of observation and several standard interpretations are 
presented in a new manner. Thereafter, the actual situations 
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on the effect of Mach number, sweepback effect and effect of 
thinness of moving disturbance volumes and disturbance surfaces 
are formulated on a new basis. The solution of the wave equa- 
tion is then described on a physical basis and the field 
quantities produced from this in a general form. The present 
unified treatment applies as well to aerodynamics and aero- 
acoustics as it does for subsonic and supersonic Mach numbers. 


33 




Fig. 1: Wave Propagation of a Disturbance Source with Smooth Motion, 
Represented in a Medium-Fixed Reference System with the Same 
Motion, Related to One Another Through the Galilei Transformation. 

Key: a. Reference System with same motion, 
b. Medium-Fixed Reference System 




Fig. 2: Kinematic Relationships in the Disturbance Pro 
pagation of Moving Disturbance Sources in a 
Homogeneous Medium at Rest with Infinite Dimen 
sions . 




Fig. 3: Universal' 'Relationships for Determination of Radiation 
Quantifies r v and 9 V of Disturbance Signals, Hitting a 

Moving Disturbance Source at Varying Mach Numbers Ma s 

at the Field Points P(r„, 9 ) at Time t_. 
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Fig. A: Production of the Spatial Correlation Functions 

L and a N as a Re sult of the Actual Longitudinal 

Sith r a?? V n^ S ? E * tension of the Emission Elements 
with all Disturbance Signals Reaching Field Point 
P Simultaneously at Time t . nt 

, _ o 


d V dl o 


= dn /dn ) 
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Fig. 5: Production of the Spatial Correlation Functions 
as a Result of the Actual Volume Expansion of 
the Emission Elements in the Case of Various 
Orientations of the Disturbance source in Space 


or to the Path Line 


o ) 

N • 











Fig. 6: Production of the Spatial Correlation Function ° 
as a Result of the Actual Shift in Emission R 
Points of Moving Source-Sink Elements of a Dipole 
with the Disturbance Signals Reaching Field Point* 
P Simultaneously at Time t 


= dr /dr ) . 
v vo 
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Fig. 7: 


Production of the Doppler Factor cr D and a* as a 


Result of the Change in Wave Passage Duration dt Q 
through a resting or moving Field Point in 
Comparison to Duration dt v of the Emission of 
Disturbance Signals. 
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Fig. 8: The Doppler Factor a D in the Propagation of Spherically 

Symmetrical Waves from Moving Disturbance Sources in a 
Homogeneous Medium at Rest. 
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Mo s -= 0 bzw.'d' v =l£ 
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Fig. 9: a Universal Representation of the Correlation 
Functions a y, Op and Op in the Propagation of 

Spherically Symmetrical Disturbance Waves from 
Moving Singularities. 



Mo s Po(U 


& 0 = 0 bzw.jr 
Oder Mas z O 


0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,6 0,9 1,0 
Ma$ sin ^ — ► 


Re P res ® nt ation of the Actual Emission 
Radii r v /r Q for Random Field Points in Space 

with Subsonic and Supersonic Mach Numbers of the 
Disturbance Sources. 





Fig. 11 : Radiation Processes of Moving Source Lines of 
Infinite Length with Various Mach Numbers and 
Various Orientations to the Path Line with the 
Disturbance Signals from the L v lines Reaching 

the Field Point P(x,0,0) Simultaneously at 
time t . 
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Ma $ =y/2' ; X= 30 
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Fig. 12: Radiation Processes from a Moving Source Line in the Case 
of Supersonic Mach Numbers and Correlated Emission Lines L 

and L with Disturbance Signals reaching Field Point P at 
limes t Ql and t Q2 and always producing Constant Field 
Quantities for x :> x_ . 
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Ma s 2:1 iy Xs = 0 



Fig. 14: Surfaces in Space with Constant Values of r* at 

v -■ 

Subsonic, Sonic and Supersonic Velocities of the 
Disturbance Sources with the Axis x of Symmetry 
Through the Field Point P. 
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60 ° 80 ° 100 ° 120 ° 



60 ° 80 ° 100 ° 120 ° 



Fig. 15: Near Field Protions of the Disturbance Potential 
on a Unit Sphere with r Q = 1 around with 

Various Mach Numbers of a Moving Source and a 
Moving Dipole. 

= vs or V» cos ^ ] ' 
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Fig. 16: Distant Field Portion of the Disturbance Potential 
in the Disturbance Waves with r v = 1 around K v at 

Various Mach Numbers of a Moving Source and a 
Moving Dipole. 

(?£ = 4> F /S or 4>p/ w^Dcos’J') . 
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Fig. 17'- Near Field and Distant Field Portions of the 
Disturbance Quantities on a Unit Sphere with 

r = 1 around £ or with r = 1 around ^ 
o o v 

with Various Mach Numbers of a Moving Source, 
(sg = s N /w v s or s* = Sp/cjJ D C osi|i ). 
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Fig. 


18: 


The Ratio of the Nonsteady to the Quasi-Steady Field Quantities 
in the Near Field and Distant Field Portions of a Moving 
Source 


( X* 


s nonst . 



. st . 
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